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Quark loop contribution to tt° — > 47 

Yu.M. Bystritskiyfl V.V. Bytev, and E.A. Kurae\|l 

Joint Institute for Nuclear Research, 14-1980 Dubna, Russia 
(Dated: February 1, 2008) 

We find the contribution of constituent quark loop mechanism to the branching 
ratio of 7r° — > 4*y process to be Bfe dr = r 7r o_ j .4 7 /r 7r o_, 2 7 ~ 5-45 • 10~ 16 for the reason- 
able choice of constituent quark mass m ~ 280 MeV. This result is in agreement with 
vector-dominance approach result obtained years ago. Thus the main contribution 



arises from QED mechanism ir° — > 7(7*) — > 7(37) including light-light scattering 
C^ - ', block with electron loop. This contribution was investigated in paper of one of us 



and gave B^ ED ~ 2.6 • 10 -11 . 



I. INTRODUCTION 

In the program of investigation of CP- violating effects one of directions is the search of the 
"forbidden" process 

Oh, 

7T° -> 3 7 . 

■ An important background to it is the allowed decay: 

U ; 7T° -f 4 7 . (2) 

A realistic estimation of the differential width of it is needed. Two types of contributions - 
hadronic and electromagnetic - as a relevant mechanisms of 7r° — > 47 must be considered. 

It was suggested in |3] (1972) and confirmed in j^J (1995) that the main contribution have 
an electromagnetic nature: 

b qed = ^ 2 6 _ 10 _u (3) 



47 r 

As for hadronic part - two quite different predictions was done: paper 



5] 



gave Bt dr ~ 10" 9 and 



paper |3j result is B^ dr ~ 10 14 . Significantly lower bound for branching B^ dr ~ 7.1 • 10 18 
was quite recently obtained in frames of chiral perturbation theory in p. 
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FIG. 1: Diagrams of hadronic contribution within Vector-Dominance model (VDM). 
This contradiction is the motivation of our investigation. 

It's known the duality property in description of strong-interaction phenomena: 

• The meson and barion (hadrons) approach; 

• Quark-gluon approach. 

In this paper we estimate (using the quark-gluon approach) the 7r° — > 47 width in the model 
with constituent quark loop. 



II. APPROACH BASED ON PCAC AND VECTOR-DOMINANCE HYPOTHESES 

Considering the process 7r° — ► 47 within the Vector-Dominance model (VDM) leads to the 
following Feynman diagram (FD) Fig. ^ The amplitude which corresponds to this FD is jjj: 



M, 



2l 



M 



9piry Quiry 



(k 2 + h + h) 2 - m 2 (k 2 + k 3 + fc 4 ) 2 - ml 



1 

x 



(k 3 + h) 2 - M 2 



x e^x e£k%elk$ ■ e aM5 (k 2 + k 3 + k 4 fejk( ■ efag, (k 2 + k 3 + k 4 f e\ fc*' + 

+ 11 terms obtained by permutations, (4) 

where {e^, k{\ - is the polarization vector and momentum of i-th photon; M — is the mass of 
pion and M 2l - is the amplitude of 7r° — > 27 decay. The coupling constants are fixed by VDM 
to be: 

_ 9ttw P K ,r\ 

y 1,717 M 3 

where A p = g p , K = 9a3 /V2, 9 2 p = 2F 2 m% gljAir = 0.51 (see F n = 94 MeV and: 

1 0.43 g 2 Das, Mathur, Okubo [8j. 

9<* = \ , r , (6) 

0.23 g 2 Oakes, Sakurai 



FIG. 2: Diagrams of hadronic contribution within constituent quark loop approach. 

Estimating the branching ratio contribution which comes from amplitude (JIJ) one gets jjj: 

[ 8.6 ■ lCr 16 , Das, Mathur, Okubo. 

B hadr < ) _ ( 7 ) 

I 7.0 • 10 16 , Oakes, Sakurai. 

In paper \^ the photons identity was not taken into account (thus destructive interferences 
were lost). And this gave the strongly overestimated result B^ dr ~ 10~ 9 . 

III. APPROACH OF CONSTITUENT (HEAVY) QUARK LOOP 

In this paper we consider the mechanism with the constituent quark loop. First we need to 
determine the pion-quark coupling constant g. To do this we write out the n° — > 27 amplitude 
within quark-loop approach assuming the pion wave function to be 

7T° = — (uu-dd). (8) 
\/2 



This gives (see Fig. a): 



M 27 = 2z</^ ^= F(z)(e 1 e 2 k 1 k 2 ), (9) 



where {e\e 2 kik 2 ) = £ fj.ua/3G1e2k" k^, a = 1/137 - is fine-structure constant, the color-charge 

factor is N 2 = 3 • ^(f) 2 — (|) 2 j — 1, m, M - constituent quark mass and neutral pion mass 

1 

correspondingly, z = M 2 /m 2 and F(z) = J — ln(l — zx{l — x)). We use constituent quark 
mass equal torn = 280 MeV (according to the analysis performed in paper |lo|). thus z ~ 0.23 
and F(z) « -0.118. 

Comparing with current- algebra (CA) result: 

m ^ a = ^r^ k ^ v ^ = 4^ w 7A ev ' (io) 
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with F n = 94 MeV, we obtain for pion-quark coupling constant g = 2.065. 

Now we start to estimate the hadronic contribution to tt° — > 47 decay within quark-loop 
approach which gives us the set of FDs presented at Fig. b. The amplitude for the first FD 
has the form: 

where A 1 = k 2 — m 2 , A 2 = (k + k 2 ) 2 — m 2 , A 3 = [k — k 3 ) 2 — m 2 , A A = [k — k 3 — /c 4 ) 2 — m 2 , 
A5 = (k+ki+k 2 ) 2 —m 2 ; k,m- are the momentum and mass of quark in the loop (m = 280 MeV). 
Sp = Sp 7 5 (/c — /c 3 — k 4 + m)e 4 (fc — fc 3 + m)e 3 (/c + m)e 2 {k + k 2 + m)ei{k + fcj + k 2 + m) . 
The color-charge factor for this case is N4 = 3 ■ ^(f) 4 — (§) 4 ) = §• Next we join the de- 
nominators of loop fermions propagators Af 

4! / * „ (12) 



Ai---A 5 J (A lXl + ■ ■ ■ + A 5 x 5 ) 


where dr = dx\ . . . dx§ S(x\ + • • • + £5 — 1), and get A1X1 + • • • + A 5 x$ = (k + b) 2 — d, where 
6=(fci + k 2 )x 5 + k 2 x 2 - k 3 x 3 - (k 3 + fc 4 )x 4 , d = m 2 - M 2 a 2 , a 2 = x 4 x 5 + x 3 a;5(ti 3 + £33) + 
t22,x 2 x 3 + x 2 x 4 (t 2 3 + £24) + £34^4(^1 + x 2 ) + £12^5(^1 + ^3), Uj = 2{k i k j )/M 2 . Performing a 
shift of loop momentum (k — > x + 6) and calculating the loop trace Sp[. . . ] we can cancel even 
powers of new integration momenta x with denominators (x 2 — d). After that we may use a 
well-known formula of integration: 

d\ 1 _ (-1)" 1 (u] 

m 2 (x 2 - d) n (n - 1) {n - 2) d n ~ 2 ' 1 ' 

To perform the integration over Feynman x's we will expand the obtained expression for am- 
plitude over z = M 2 /m 2 ~ 0.23, i.e. use formula: 

dr n = m~ 2n ^1 + na 2 z + ^n(n + l)(a 2 z) 2 + ^n{n + l)(n + 2){a 2 zf + . . . J . (14) 

This makes the integration over x's trivial: 

dT T n lT n 2 n z n, _ ^\ U 2 \ U 3 \ n 4 ! 

dr Xl x 2 x 3 x 4 ~ {ni+n2 + n3 + n4 + 4)] - W 
Thus we obtain the amplitude for first FD: 

Mff = -ga 2 ^= • /, I = - (B + zB 1 + z 2 B 2 + z 3 B 3 + ...) , (16) 
V 2 Tn. 
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where Bi - some rather complicated but known functions of and tij. The total amplitude 
M 47 is the sum of 24 terms obtained from (fTBj) by permutations of final photon legs (i.e. pairs 

M 47 =^M«. (17) 

perm 

Performing this summation we see that: 

J2 Bo = J2 Bi = E B * = °> J2 B ^°- ( 18 ) 

perm perm perm perm 

This can be obtained by applying Schouten identity: 

{pxP2pm)Qv, = (vP2P3Pd(Qpi) + ipiwm){Qp2) + (pmw^iQps) + {pxP2Pzv){Qp±), (19) 

to antisymmetric Levi-Civita tensors (aft'yS) in functions Bi. Here we should notice that total 
amplitude M 47 fulfills the requirements of gauge invariance and Bose symmetry: 

K Mp up , = K Mp up , = K MX, = K MX, = 0, 

MXaih, h, h, h) = MXAh, h, k 3 , k 4 ) = ... . 

After squaring the amplitude M 47 and summing over photons polarizations {AiA 2 A 3 A 4 } we 
have: 



E |M 47 ^(- 9a ^) 2 ^iip, with p =j:s 2 , 

A1A2A3A4 perm 



(20) 



where: 



S 2 = ~ti 2 {8t 13 (t 2 u (t 12 t 13 + t 23 (-2t u + 3t 23 )) + 

+ (3tl2^13 + 4^23 ( — 2t 44 + £23)) ^24 + 6^13^24) + 4^13^14 (^14 — ^23) ^24^34 — 

- 32£ 13 44 + (3^ 2 - 40t 13 t 2 4) 4 + 13ti 2 4) • (21) 
The 7T° — > 47 decay width dT^? has a form: 

dr " r= (2^ml £ |M ^ 2 *' (22) 

v 7 ^ A1A2A3A4 

where the phase volume <i0 4 is considered in Appendix [XJ To evaluate the hadronic contribution 
to total decay width T 47 dr we perform the numerical phase-space integration and finally obtain: 



47 

T2 



M N 2 f 

= 2i%6 g2a '^t zlJ) Where J = J d ^ P ^ °- 097 ' ( 23 ) 

where the dimensionless phase volume <i0 4 is defined by relation <i0 4 = 1 yM 4 (i0 4 . Thus our 
result for hadronic contribution to total decay width is T^ dr « 4 ■ 10~ 15 eV. 
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IV. CONCLUSION 

Using the result listed above (see (|23p) we estimate the hadronic contribution to branching 
of 7T° — > 47 decay as: 

j^hadr 

B hadr = ^ . 1Q -16 ^ 



This value has the same order of magnitude as one obtained in [3j . It also is in a good agreement 
with the recent result llll] obtained within the meson and barion approach. But it contradicts 

n 

the result of paper 6] obtained by using chiral perturbation theory. 

The result obtained in MJ B%f r < 10" 9 is in strong contradiction to our estimate and QED 
. As for result of paper [5J] it presumably needs a serious revision as well as amplitudes 
interferences was not taken into account. 

We should notice that our consideration were performed within QED. As for QCD corrections 
they may be parameterized in the variation of constituent quark mass m. For instance if 
uncertainty of m is about 5 MeV it can result in SB^/B^ ~ 30%. So these corrections can be 
essential. 

Finally we may conclude that the main contribution to 7r° — > 47 decay has the QED nature, 
i.e. mechanism with light-light scattering block with electron loop: tt° — > 7(7*) — > 7(37). This 
contribution was investigated in paper of one of us [4] and appeared to be equal to 

B^ ED w 2.6- KT 11 . (25) 



APPENDIX A: PHASE VOLUME 

The phase volume dcj)^ for decay tt° — ► 47 is: 

LiO\ IUJ2 2^3 ZUJ4 

Being expressed in terms of energy fractions of photons jji = uJi/M (which satisfy the identity 
Hi + 2/2 + U3 + Hi = 1) and cosines of photons momenta orientations CV,- = cos(k{kj) this phase 
volume reads to be: 

2 1/2 2tt 1 

J 'dfa ... = y M 4 J ' d Vl dy 2 J dfo J dC l2 dC 13 ^^ . . . , (A2) 



-1 
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where A = l — y 1 (l — C r3 ) — y 2 (l — C 23 ). The region of integration is determined by conditions: 
< 2/1,2 < ^ - 1 < °iJ < 

< 2/3 = ^(1 - 2(y! + |/ 2 ) + 2y m (l - C 12 )) < ^ (A3) 
< 2/4 = \Jvl + yl + yl + 2ymC 12 + 2y 1 y 3 C 13 + 2y 2 y 3 C 23 < ~, 

while C 23 = C12C13 + Si 2 Si 3 cos(0 3 ), and = sin(fcjfcj). The kinematical variables £jj can be 
expressed as: 

£l2 = 2 2/12/2(1 — C12), ^23 = 2 2/2 2/3 (1 - C 23 ), 

t 13 = 2 2/x 2/3 (1 - C13), £24 = 22/2(2/2 + 2/4 + 2/1^12 + 2/3^23), ( A4 ) 

*14 = 2 2/i (2/1 + 2/4 + 2/2C12 + 2/3C13), ^34 = 2 2/3 (2/3 + 2/4 + 2/1C13 + 2/2^23)- 
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